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Abstract. Under general conditions, the equation g(x,y) = implic- 
itly defines y locally as a function of x. In this article, we express divided 
differences of y in terms of bivariate divided differences of g, generalizing 
a recent result on divided differences of inverse functions. 



1. Introduction 

Divided differences can be viewed as a discrete analogue of derivatives 
and are commonly used in approximation theory, see [l] for a survey. 

Recently, the second author and Lyche established two univariate chain 
rules for divided differences [2], both of which can be viewed as analogous 
to Faa di Bruno's formula for differentiating composite functions (6|j7|. One 
of these formulas was simultaneously discovered by Wang and Xu fSj. In a 
follow-up preprint, the other chain rule was generalized to the composition of 
vector- valued functions of several variables [4] , yielding a formula analogous 
to a multivariate version of Faa di Bruno's formula fs]. 

In [3 1, the univariate chain rule was applied to find a formula for divided 
differences of the inverse of a function. In Theorem [T| the Main Theorem 
of this paper, we use the multivariate chain rule to prove a similar formula 



for divided differences of implicitly defined functions. Equation 16 shows 
that the formula for divided differences of inverse functions in ^ follows as 
a special case. 

More precisely, let y be a function that is defined implicitly by a function 
51 : — )• M via g[x,y{x)) = and ^[x,y{x)^ ^ 0, for every x in an open 
interval i7 C M. Then the Main Theorem states that for any 

xo, . . . , G C/, yo := ^(a^o), • • • , := y(xn) G y{U) 

we can express [xq, ■ ■ • , x„]y as a sum of terms involving the divided differ- 
ences [xio , . . . , Xj, ; j/i, , . . . yi^\g, with < < n < • • • < ir < n- 

In Section [2| we define these divided differences and explain our notation. 
In Section [3j we apply the multivariate chain rule to derive a formula that 
recursively expresses divided differences of y in terms of divided differences 
of g and lower order divided differences of y. Finally, in Section [4j we solve 
this recursive formula to obtain a formula that expresses divided differences 
of y solely in terms of divided differences of g. We end the section with 
applying the Main Theorem in some special cases. 
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2. Divided Differences 

Let [xq, . . . , Xn]f denote the divided difference of a function / : (a, 6) — t- M 
at the points xq, . . . ,Xn, with a < xq < ■ ■ ■ < Xn < b. If all inequalities are 
strict, this notion is recursively defined by [xq]/ := /(xq) and 

r 1. [xi,...,Xn]f -[xo,...,Xn-l]f ^„ 

[Xq,. . . ,Xn\f = lfn>0. 

Xfi Xq 

If some of the {xj} coincide, we define [xq, • • • , x^]/ as the limit of this for- 
mula when the distances between these {xj} become arbitrary small, pro- 
vided / is sufficiently smooth there. In particular, when xq = • • • = x„, one 
can show that [xq, • • • ,x„]/ = f^'^\xQ)/n\ . For given io, . . . ,ik satisfying 
^0 ^ ^1 ^ ■ ■ ■ ^ ifc) we shall sometimes shorten notation to 

(1) [ioh ■ --iklf ■■= [xio,Xi,, . . .XiJ/. 

The above definitions generalize to bivariate divided differences as follows. 
Let / : [/ — 7- M be defined on some 2-dimensional interval 

U = (ai,6i) X (02,62) C M^. 

Suppose we are given m,n > and points xq, . . . , Xm £ (oi, 61) satisfying 
Xq < ■■■ < Xm and yo,...,ym G (02, 62) satisfying yo < ■ ■ ■ < Vm- The 
Cartesian product 

{xo, . . ■,Xm} X {yo, ...,yn} 

defines a rectangular grid of points in U. The (bivariate) divided difference 
of / at this grid, denoted by 

(2) [xo,...,Xm;yo,---,yn\f, 

can be defined recursively as follows. If m = n = 0, the grid consists of only 
one point (xo,yo)) and we define [xo;yo]/ := f{xo-,yo) as the value of / at 
this point. In case m > 0, we can define ^ as 

[Xi, . . . , Xn^; j/o, . . . , yn]f - [xq, . . . , Xpt-i; j/Q, . . . , yn]f 

1 

Xm Xq 

or if n > 0, as 

[xq, . . . , x^; • • ■ , yn]f - [xq, . . . , Xn,; j/0, • • • , 
Vn-yo 

If both m > and n > the divided difference ^ is uniquely defined by 
either recursion formula. 

As for univariate divided differences, we can let some of the points coalesce 
by taking limits, as long as / is sufficiently smooth. In particular when 
Xo = • • • = Xm and yo = • • • = yn, this legitimates the notation 

Qm+n J 

[xo, ...,Xm;yo,-- • ,yn]/ := — r-, ^ ma v i^o,yo)- 

mini ox"^oy^ 

Similarly to Equation [T| we shall more often than not shorten the notation 
for bivariate divided differences to 

(3) [io^i ■■■is; joji ■■■jt]f-= [xio , Xii , . . . , Xi^ ; yj^ , , . . . , 
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3. A Recursive Formula for Implicit Functions 

Let y be a function implicitly defined by g[x,y{x)) = as in Section ol 
The first step in expressing divided differences of y in terms of those oig 
is to express those of g in terms of those of y. This link is provided by a 

special case of the the multivariate chain rule of | Let M A ^ M 
be a composition of sufficiently smooth functions f = ((/>, ip) and g. In this 
case, the formula of for n > 1 is 

n k 

(4) [xo,xi,...,2;„](5of) = ^ ^ ^ 

fe=l 0=io<ii<---<jfc=n s=0 

[<P{xio),<P{xi^), ■ ■ . ,(l){xiy,%l){xi^),i){xi^^^), . . .,il}{xij\g 

s k 
1=1 l=s+l 

Now we choose f to be the graph of a function y, i.e., f : x i— )• ((/)(x), = 
(x,y(x)). Then the divided differences of cj) of order greater than one are 
zero, implying that the summand is zero unless (io, h, ■ ■ ■ ,is) = (0, 1, ... , s); 
below this condition is realized by restricting the third sum in Equation |4] 
to integers s that satisfy s = is — io- Since additionally divided differences 
of (p of order one are one, we obtain 

n k 

(5) [xo,xi,...,x„]5(-,y(-)) = ^ ^ J] 

k=l 0=io<ji<---<ife=n s=0 

s=is—io 

k 

[xq, xi, . . . , Xs', yig, yig^i, ■ ■ ■ , yif.]g [xjj_j^, Xjj_j_|_i, . . . 

l=s+l 

where yj := y{xj) for j = 0, 1, . . . , n. For example, when n = 1 this formula 
becomes 

[xo, xi]g{-,y{-)) = [xo; yo, yi]g [xq, xi]y + [xq, xi; yi]g, 
and when n = 2, 

[xo,xi,X2]g{-,y{-)) = [xo;yo,y2]g[xo,xi,X2]y 

+ [xo;yo,yi,y2]g [xo,xi]y [xi,x2]y 
+ [xo,xi;yi,y2]g [xi,X2]y 
+ [xo,xi,X2;y2]9- 



In case y is implicitly defined by g(x,y(x)) = 0, the left hand side of 
Equation [5] is zero. In the case n = 1, therefore, we see that 

where we now used the shorthand notation from Equations [T] and [3j For 
n > 2, the highest order divided difference of y present in the right hand 
side of Equation [5] appears in the term [0; On]g [01 • • • n]y. Moving this term 
to the left hand side and dividing by —[0;On]g, one finds a formula that 
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expresses [01 • • • n]y recursively in terms of lower order divided differences of 
y and divided differences of g, 

(7) [oi...%=-x: E E 

k=2 0=io<---<*fc=n s=0 

s=is—io 



k 

- ■ ■ s]isis+i- ■ ■ik]9 -TT r. f. , ,s -1 

^ n [n-i{^i-i + i)---^i]y. 

^ ' -'^ l=s+l 

We shall now simplify Equation [7| By Equation [6j the first order divided 
differences of y appearing in the product of Equation [7] can be expressed 
as quotients of divided differences oi g. To separate, for every sequence 
(io, ii, . . . , ifc) appearing in Equation [Tj the divided differences of g from 
those of y, we define an expression involving only divided differences of g, 

k k / 

r. . . [io---is;is---ik]9 tt / [ii-iii;ii]g 

fr^Q h;ioik]g ^^^^^ \ [ii-i;ii-iii]g 

s=is—io ij— j;_j=l 

Note that if a sequence {iq, . . . ,ik) starts with precisely s consecutive inte- 
gers, the expression {iq ■ ■ ■ ik}g will comprise s terms. For instance, 



{013}g 



[0;03]g [2;23W 
[0;013]g[01;l]g [01; 13]ff 
[0;03]ff [0;01]5 [0;03]5' 



lOm [0;012]g[01;l]g[12;2]g [01; 12]g [12; 2]g [012; 2]g 

^ [0;02]5 [0;01]5[1;12]5 [0; 02]<7 [1; 12]^ [0;02]5- 

The remaining divided differences [ii-i ■ ■ ■ ii]y in the product of Equation 
[7] are those with ii — > 2, and each of these comes after any s satisfying 
s = is — iQ- We might therefore as well start the product of these remaining 
divided differences at Z = 1 instead of at I = s + 1, which has the advantage 
of making it independent of s. Equation [7] can thus be rewritten as 



[0---ri]y = E {io---ik}g n 



k=2 0=io< — <ik=n 1=1 

For n = 2,3,4 this expression amounts to 

(9) [012]y={012}5, 

(10) [0123]y = {0123}<7 + {023}^ [012]?/ + {013}5 [123]y, 

(11) [01234]y = {01234}^ + {0134}5 [123]y + {034}^ [0123]y 

+ {0124}5 [234]y + {0234}^ [012]y + {014}g [1234]y 
+ {024}5 [012]y [234]y. 
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4. A Formula for Divided Differences of Implicit Functions 

In this section we shall solve the recursive formula from Equation [7[j 
Repeatedly applying Equation [tJ] to itself yields 

(12) [012]y={012}5, 

(13) [0123]y = {0123}5 + {023}^ {^I2}g + {013}g {123}<7, 

(14) [01234]y = {01234}g + {QlM}g{l2^}g + {QM}g {013}^ {123}5 

+ {034}5 {0123}5 + {034}5 {023}^ {Ql2}g + {0124}5 {234}^ 
+ {0234}5 {012}5 + {014}5 {134}5 {123}5 + {014}^ {1234}5 
+ {014}5 {124}5r {234}<7 + {024}^ {012}^ {234}^. 

Examining these examples, one finds that each term in the right hand sides 
of the above formulas corresponds to a partition of a convex polygon in a 
manner we shall now make precise. 

With a sequence of labels 0, 1, . . . , n we associate the ordered vertices of a 
convex polygon. A partition of a convex polygon is the result of connecting 
any pairs of nonadjacent vertices with straight lines, none of which intersect. 
We refer to these straight lines as the inner edges of the partition. We 
denote the set of all such partitions of a polygon with vertices 0, 1, . . . ,n 
by 7^(0, 1, . . . , n). Every partition vr G 7^(0, 1, . . . , n) is described by its set 
F[-k) of (oriented) faces. Each face / G -F'(7r) is defined by some increasing 
sequence of vertices io, ii, • • . , ifc of the polygon, i.e., / = (zq, ^i, • • • , ifc)- We 
denote the set of edges in tt by E{tt). 

Let y be a function implicitly defined by g[x, y(x)) = and (xq, yo), . . . , 
{xn^Un) be as in Section [TJ Equations 12 -14| suggest the following theorem. 



Theorem 1 (Main Theorem). For y and g defined as above and sufficiently 
smooth and for n >2, 

(15) [0---n]y= Yl n {vo---^r}g, 

7reP(0,...,n) (t)o,...,Dr)GF(7r) 

where {vq ■ ■ ■ Vr}g is defined by Equation^ 

Before we proceed with the proof of this theorem, we make some re- 
marks. For n = 2,3,4: this theorem reduces to the statements of Equations 
[T2HHI To prove Theorem [TJ our plan is to use Equation [tJJ recursively to 
express [01 ••• n]y solely in terms of divided differences of g. We have found 
it helpful to assign some visual meaning to Equation [TjJ Every sequence 
i = {iQ,ii, . . . ,ik) that appears in Equation [tJJ induces a partition vTi S 
V{0, 1, . . . , n) whose set of faces comprises an inner face {io,ii, . . . ,ik) and 
outer faces {ij,ij + 1, . . . , ij+i) for every j = 0, . . . , k — 1 with ij+i — ij > 2. 
We denote by "Pi the set of all partitions of the disjoint union of these outer 
faces. An example of such a sequence i, together with its inner face, outer 
faces, and partition set Vi is given in Figure [TJ 

We shall now associate divided differences to these geometric objects. To 
each outer face {ij,ij + l, . . . , ij+i) we associate the divided difference [ij{ij + 
1) • • • ij+i]y, and to each inner face {io, ii, . . . , Zfc) we associate the expression 
{^0 ■ ■ ■ 'ik}g- For any sequence i that appears in the sum of Equation!? J the 
corresponding inner face therefore represents that part of Equation [7| that 
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3 



2=1 




0=1(1 7=^5 



Figure 1. For n = 7, the sequence i = (0, 1, 2, 4, 5, 7) gives 
rise to the two outer faces (2,3,4) and (5,6,7), which are 
drawn shaded in the figure. The set Vi contains in this case 
just 1x1 = 1 partition, namely the union of the unique 
partitions {(2, 3, 4)} and {(5, 6, 7)} of the outer faces. 

can be written solely in terms of divided differences of g, while the outer 
faces represent the part that is still expressed as a divided diff^erence of y. 

Repeatedly applying Equation [tJ yields a recursion tree, in which each 
node represents a product of divided difference expressions associated to 
inner and outer faces. These recursion trees are depicted in Figure [2] for 
n = 2, 3, 4. Equation [tJ roughly states that the expression of any nonleaf 
vertex is equal to the sum of the expressions of its descendants. 

Proof of the Main Theorem. This theorem is a generalization of Theorem 1 
in |3], and the proofs are analogous. We prove the formula by induction on 
the order n of the divided difference of y. 

By the above discussion, the formula holds for n = 2, 3, 4. For n > 5, 
assume the formula holds for all smaller n. Consider the recursive formula 
from Equation [tJ} For every sequence i that appears in this equation, the 
corresponding outer faces have fewer vertices than the full polygon. By 
the induction hypothesis, we can therefore replace each divided difi^erence 
[ir ■ ■ ii-\-i]y appearing in the product of Equation [tJ by an expression involv- 
ing only divided differences of g. 

As before, let Vi denote the set of all partitions of the disjoint union of the 
outer faces induced by i. Then, by the induction hypothesis, the product in 
Equation [7| is equal to 



For a given inner face i, the set Vi can be identified with {vr E V{0, . . . ,n) : 
i G F{7r)} by the bijection F{7r) i— )• -F(7r) U {i}. Substituting the above 
expression into Equation (tI] then yields 





TT&Vi {vo,...,Vr)eF{TT) 





i=(io,...,ifc) 



inner faces 



TTSPi (vo,...,Vr}l^F{Tr) 






inner faces •7rGP(0,...,n) {vo,...,Vr)£F{n) 
i={io,...,ik) ieF(7r) 




□ 



7reP(0,...,n) (vo,...,Vr)eF(Tr) 
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{0123} 







Figure 2. For n = 2,3,4, the figure depicts the recursion 
trees obtained by repeatedly applying Equation [7fl The top 
levels of these recursion trees correspond to Equations [9 - 11 



Intuitively, this proof can be expressed in terms of the recursion tree as 
follows. As remarked in the previous section, Equation [7|] states that the 
expression of any nonleaf vertex is equal to the sum of the expressions of 
its descendants. By induction, the expression [01 • • • n]y of the root vertex 
is therefore equal to the sum of the expressions of the leaves, which, by 
construction, correspond to partitions of the full polygon. 
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Example 1. Let us apply Theorem [T] to fin d a simple expression for divided 
differences of the function y{x) = ^/l — defined on the interval (—1,1). 
This function is implicitly defined by the polynomial g{x, y) = x^+y^ — 1 = 0. 
For any knots Xa,Xb,Xc,Xd satisfying — 1 < Xa < < Xc < < 1 and 
corresponding function values ya,yb,yc,yd, one finds 

[xa, Xb] ydg = Xa + Xb, [xa, Xf,, Xc, yd]g = 1, 



[xa;yb,yc\g = yt + yc, 



[xa;yb,yc,yd]g = l, 



and all other divided differences of g of nonzero order are zero. In particular, 
every divided difference of g of total order at least three is zero, which means 



that the sum in Equation 15 will only be over triangulations (i.e., partitions 
in which all faces are triangles). For a polygon with vertices 0,1,..., n. 
Exercise 6.19a of 12 states that the number of such triangulations is given 



by the Catalan number 



C{n-1) 



2n 

n - 



1 



Consider, for a given triangulation vr G 7^(0, 1, . . . , n), a face (a, 6, c) S 



-F(7r) from the product in Equation 15 As any divided difference of the 
form [xa, Xf,; yh, yc]g is zero for this g, Equation [s] expresses {abc}g as a sum 
of at most two terms. There are four cases. 



{abcjg 



-1 



ya + yc 
1 

ya + yc 
ya +^yc 
ya + yc 



1 + 



ya + yb 

Xa ~l~ Xb 

ya + yb 
yb + yc 



yb + yc 



a,b,c consecutive; 
only a, b consecutive; 
only b, c consecutive; 
otherwise. 



For example, when n = 3, our convex polygon is a quadrilateral, which 
admits C(3 — 1) = 2 triangulations vri and tt2 with sets of faces 

F(7ri) = {(0, 1, 2), (0, 2, 3)}, F{7T2) = {(0, 1, 3), (1, 2, 3)}. 

One finds 

[xo, XI, X2, X3]Vl-x2 = {012}^ {023}g + {013}^ {123}(7 = 



-1 



(yo + y3)(yo + y2) 



{yo + y3){yi + y3) 



1 + 



1 + 



Xo + Xi Xi + X2 

yo + yi yi + 2/2 

Xi + X2 X2 + X3 

yi + y2 2/2 + ys 



X2 + X3 

y2 + 2/3 

Xo + Xi 

yo + yi' 



Example 2. Next we show that Theorem [T] is a generalization of Theorem 
1 of J3j, which gives a similar formula for inverse functions. To see this, we 
apply Theorem [T] to a function y implicitly defined by a function g{x,y) = 
x — h{y). Referring to Equationlsl we need to compute [iq - ■ ■ is] is - ■ ■ ik]g for 
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this choice of g and various indices io, . . . ,ik and s £ {0, . . . ,k}. Applying 
the recursive definition of bivariate divided differences, one obtains 



^10 



if s = 0, s = k; 
[io - ■ ■is;'is - ■ - iklx = { 1 ifs = l,s = k; 

otherwise, 



[is---ik\h ifs = 0; 
otherwise. 



[io - ■ ■is;is- ■ - ikMy) -- 
Consider a face / = (vq, . . . , Vr) of a given partition vr G V{0, . . . ,n) in 



Equation 15 Since r > 2, the divided difference [vo ■ ■ ■ Vs;vs ■ ■ ■ Vr] {x — h{y)) 
is zero for s > 1. Using this, Equation [s] expresses {vo ■ ■ ■ Vr}g as a single 



term 



r T [vo;vo---Vr]g Vt f h^ivi;vi]9 
{vo---Vr}g = , — M 



[vo;voVr]g fj^ \ [vi-i;vi^ivi]g 

[vo ■ --Vrjh 



n ^ . 

[voVr]h ^1 [vi^ivi]h 

Vl-Vl_i = l 

Taking the product over all faces in the partition vr, the denominators of the 
factors in the above equation correspond to the edges of the partition, while 
the numerators correspond to the faces of the partition. As there is a minus 
sign for each face in the partition, we arrive at the formula 

[voVl---Vr]h 

(16) [01... n]y= (-1)#^W ^"-"'-'"^^^-^ , 

{vo,vi)£E{n) 

which appears as Equation 11 in |3]. 

Note that the inverse of the algebraic function y = \J\ — in Example [l] 



is again an algebraic function. Equation 16 would therefore not have been 
of much help to find a simple expression for divided differences of y. In 
fact, Example [T] can be thought of as one of the simplest examples for which 
Theorem [T] improves on Equation 16 as it concerns a polynomial g with 
bidegree as low as (2,2). 

Example 3. In this example we shall derive a quotient rule for divided 
differences. That is, we shall find a formula that expresses divided differences 
of the quotient y = P{x)/Q{x) in terms of divided differences of P and of 
Q. Let g{x,y) = Q{x)y — P{x). Then, in Equation|8j 

yis[io- ■ -isjQ - [io- ■ ■is]P if s = k; 
(17) [io- - -is;is- - -ik]9 = { [io---is]Q if s = k - 1-, 

otherwise. 



In Equation 15 , therefore, the only partitions with a nonzero contribution 
are those whose faces have all their vertices consecutive, except possibly 
the final one. In particular, the inner face with vertices = io < - - - < 
ik = n should either be the full polygon, or should have a unique inner 
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edge {ik-i,n). By induction, it follows that the partitions with a nonzero 
contribution to Equation [15] are precisely those for which all inner edges end 
at n. These partitions correspond to subsets / C {1, 2, . . . , n — 2}, including 
the empty set, by associating with any such I the partition with inner edges 
{(i,n) : i E /}. Equation 15 becomes 

(18) [0---n]^ ={()■■ ■n}g+ 

n—2 r k 

X] {r---n}g JJ{ij_i • • • i^njc/, 

r=l fc=l 0=io<ii<---<iu=r j=l 

where the dots represent consecutive nodes and an empty product is under- 
stood to be one. A long but straightforward calculation involving Equations 
[8j[l7j and [18] yields 

r=l ^'^ k=l 0=io<ii<--<ik=r j=l ^h-i 

where Qi := Q{xi) for i = 0, . . . ,n. Alternatively, this equation can be found 
by applying a univariate chain rule to the composition x i— )• Q{x) i— )• 1/Q{x), 
as described in Section 4 of [2]. 

Finally, we note that taking the limit xq, • • • , x„ — )• x in Equations [6] [8] 
[121 and IT3] yields 



y'{x) 



910 
901 '' 



"( ^ — _ ^ _i_ n 9ll9l0 _ 9029io 

501 9oi 9oi 
III, ^ _ 530 921910 o5205ii ^ 920910902 o5i25io 

y \X j — -rOrj+Or, o o '^'^ 

901 9oi 9oi 9oi 9oi 

f. 9ll9W , glog03 , f^ giigiog02 ^ 9io9o2 
— o — 3 1 J h y J 6 — g — , 

^01 9oi ^01 501 



where we introduced the shorthand 

Qs+tg 



9st 



dx^dy* 



{x,y{x)). 



These formulas agree with the examples given in 10 
with a formula stated as Equation 7 in [ol. 



11, Page 153] and 
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